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We find in theories with spontaneous P and CP violation that symmetries needed to set the tree 
level strong CP phase to zero can also set all non-zero tree level CP violating phases to the maximal 
value 7r/2 in the symmetry basis simultaneously explaining the smallness of 9 and the largeness of 
the CKM CP violating phase. In these models we find the one loop lower bound d > lO"^'^ relevant 
for early discovery of neutron EDM d„ > 10~^^ ecm. The lower bound relaxes to ^ > 10"^"^ or 
dn > 10~^^ ecm for the case where the CP phases are non-maximal. Interestingly the spontaneous 
CP phase appears in the quark sector, not the Higgs sector, and is enabled by a heavy left-right 
symmetric vector like quark family with mass M. These results do not vanish in the decoupling limit 



of M„+ > M 



oo (where M„+ is the mass of heavy Higgses at the parity breaking scale) and 

the age-old expectation that laws of nature (or its Lagrangian) are parity and matter-antimatter 
symmetric may be testable by the above predictions and EDM experiments, even if new physics 
occurs only at see-saw, GUT or Planck scales. There is also a region in parameter space with 
M„+ < M where the above bounds are dampened by the factor (A/„+/M)^. By using flavour 

^2 2 

symmetries and texture arguments we also make predictions for the CKM phase that arises from 
the maximal phase on diagonalization to the physical basis. There are no axions predicted in this 
model. 



I. INTRODUCTION 

With the discovery of Higgs or Higgs-Uke boson the 
only standard model parameter that remains to be deter- 
mined is the value of the strong CP phase 9 that requires 
the violation of both left-right (or parity P) and matter- 
antimatter (CP) symmetries. Over the past 6 decades 
searches for the neutron electric dipole moment (EDM) 
which is also both P and CP odd have provided the ex- 
perimental bound dn ^ 1.9 x 10"^'^ ecm at 90% C.L. Ij. 
With d„ - 2 X lO-i^e* [2] this translates to 6* < 10""^°. 
Several efforts [3jare currently underway to improve the 
sensitivity of these experiments by two orders of magni- 
tude that can potentially probe d„ down to 10~^^ecm or 
6 ^ 10~^^. Storage ring experiments [4] being planned 
to search for EDMs of the proton and deutron can fur- 
ther this frontier to an equivalent of d„ ~ lO^^^ecm or 
^-10-13. 

Since both P and CP are broken in nature we would 
naively expect ^ ~ 7 where 7 is the CP violating CKM 
phase. However experimentally 6 << j ^ 69° which 
hints at a hidden symmetry, and this large inequality is 
the well-known strong CP puzzle. 

The most popular solution to the strong CP problem 
is the Peccei-Quinn (PQ) symmetry [5] that dynamically 
sets 9 to zero. The neutron EDM induced radiatively in 
the standard model in this case is dn ~ 10~^^ecm and 
is too small to be observed in the ongoing experiments. 
There can be beyond the standard model contributions to 
dn due to new physics such as supersmmetry but these 
decrease quadratically as the scale of new physics be- 
comes large. For example with supersymmetry as new 
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physics, dn ^ 10 '^'^{TeV/Msusv)'^ ecm [5] and becomes 
smaller than the sensitivity of ongoing experiments for 
MsusY > (10-100) TeV. 

In this work instead of using the PQ symmetry we 
follow a different line of approach to solve the strong CP 
problem and show that the model presented in this paper 
(along with that in reference 6 ) can lead to a sizable 
dn that is discoverable in the ongoing EDM experiments 
even if all new physics occurs only at very high scales 
such as GUT or Planck scales. 

Since violation of both P and CP is fundamental to 
the existence of a non-zero d„ the approach different from 
the continuous PQ symmetry is to impose either of these 
discrete symmetries to set 9 to zero [BHH] ■ The challenge 
then is that P/CP must break without spontaneously 
generating a strong CP phase at the tree-level so that the 
experimental constraint ^ << 7 is respected. In the past 
this has typically required along with vectorlike quarks 
further symmetries that generate the Nelson-Barr form 
of mass matrices as in [7j or mirror parity defined so 
that it takes known quarks and leptons to 3 additional 
generations of mirror quark and lepton families as in [5] . 

However more recently in [6] we showed that with no 
other symmetries imposed and no multiplication of exist- 
ing families by three generations of mirror families, just P 
and CP are sufficient to solve the strong CP problem in 
the left-right symmetric model with the addition of a full 
vectorlike quark family. A crucial aspect of this model is 
that terms and vacuum expectation values (VEVs) that 
violate P conserve CP and those that violate CP conserve 
P and therefore the strong CP phase is not generated at 
the tree level as it is protected by either P or CP for every 
term. Our work in [^ is the first solution of the strong 
CP problem in the left-right symmetric model where P 
(and not PQ symmetry)is used to set the tree-level 9 to 
zero. (See also [Hj.) 



Moreover 9 generated radiatively at the one-loop level 
in this model does not diminish as the scale of P and 
CP breaking (that is mass scale M^+ of Higgs sector 
that breaks P and the mass M of vectorlike quarks that 
are needed to break CP) go to infinity, since CP phases 
generated in collusion with Yukawa terms do not respect 
the decoupling theorem, dn generated in this model can 
be much greater than the naive standard model expecta- 
tion of 10~^^ec'm even if there is no new physics at TeV 
or lOOO's of TeV scales and therefore should be of interest 
to the ongoing neutron EDM experiments. 

We first generalize the model of reference [B] to allow 
for spontaneous CP violation (instead of softly through 
dimension 3 fermion mass terms as in that work) so that 
we can consider both the case of spontaneous (which 
could also be more predictive) as well as soft CP break- 
ing. 

To our surprise we find that thus solving the strong 
CP problem not only determines the strong CP phase, 
but also determines the CP phases in the quark mass 
matrices to be maximal (ie) tt/2 in the symmetry basis, 
consistent with the high value of the CKM CP violating 
phase 7 ~ 68° obtained on diagonalizing to the physical 
mass basis. While several works such as [IOl - fT2] study- 
ing the texture of the quark mass matrices have in the 
past suggested that CP violation could be maximal (ei- 
ther in the mass matrix or in the Jarlskog invariant), the 
phase 7r/2 is usually put in by hand and not obtained by 
symmetries. 

Spontaneous breaking of CP is an attractive idea as 
it makes CP phases calculable. In practice however the 
phases are determined by minimizing the Higgs poten- 
tial that has several additional parameters due to which 
we lose predictivity. What we show is that once we have 
both P and CP symmetries imposed and solve the strong 
CP problem, the phase of a CP violating VEV is deter- 
mined by its P transformation properties independent of 
the Higgs parameters and is maximal in the quark mass 
matrix (that also has vectorlike quarks) in the symmetry 
basis. See sections [III IHI A|IH B] and [lirc| 

It turns out that presence of a maximal CP phase 
may be experimentally verifiable by the neutron EDM 
searches. The lower bound on the strong CP phase gen- 
erated depends on rotations needed to go from the sym- 
metry basis to the physical mass basis. If the CP phase 
generated in the symmetry basis is 7r/2 (as opposed to 
an arbitrary number that could be chosen to be the ob- 
served CKM phase 7 for the purposes of calculating a 
lower bound) then some amount of additional rotation 
that mixes the real and purely imaginary terms of the 
mass matrices is needed to obtain the observed 7 ~ 68°. 
In a large region of parameter space, this generates a 
higher lower bound 9 > 10~^^ (or d„ > 10~^^ec77i) and 
the lower bound would be a couple of orders of magni- 
tude less {9 > 10-13 or d„ > IQ-^^ecm) had the CP 
phase been arbitrary and not maximal. 

The idea of the strong CP problem hinting at a hidden 
P and CP symmetry that are spontaneously or softly 



broken as presented in this work (and in [B]) is verifiable 
by finding dn greater than either of the above two lower 
bounds for the case where the heavy higgses are heavier 
than the vector like quark masses (ie) M^+ > M . While 

there is a suppression factor of {M^+ /MY if M^+ < M 
and the lower bounds are then reduced and depend on 
this factor as well. See sections llV Al and IIV Bl 

In section |V] we show how texture considerations and 
flavour symmetries can be used to make predictions for 
the CKM phase that is generated from the maximal 
phase on diagonalization from the symmetry basis to the 
physical basis. 

Finally in sections |VI| and |VII| we present some com- 
ments and concluding remarks. 



II. P AND CP PROPERTIES OF e = VACUUM 
MAXIMAL CP VIOLATION 

We assume P and CP are good symmetries that are 
both broken by VEVs of a set of Higgs fields. In gen- 
eral the vacuum will be made of several Higgs fields with 
VEVs that conserve at the tree level: 

1. both P and CP 

2. CP (but not P) 

3. P (but not CP) 

4. neither P nor CP. 

Since 9 is both P and CP odd we would expect a strong 
CP phase to be generated at the tree level by VEVs that 
conserve neither P nor CP. Now if there is a solution to 
the strong CP problem that sets the tree level strong 
CP phase to zero, we would expect that there are no 
VEVs of the fourth category. That this is in fact the 
case is proved in the next section for a class of models 
by examining the minimum of the Higgs potential. This 
means in a class of models, we can visualize the strong 
CP solving vacuum to be made of several states, with 
each state being either P even or CP even or both. This 
property of the vacuum can help determine the phase of 
CP violation that is generated at the tree level as we now 
show. 

Consider the neutral component 0° of a Higgs fields 
that picks up a VEV v. Let us say under CP (j)° — > (j)°* 
so that any non-real v breaks CP. Likewise we choose the 
CP properties of quarks and leptons so that all coupling 
constants including the Yukawa coupling are real due to 
CP. 

However as discussed, v conserves P if it violates CP 
so as not to generate a tree-level strong CP phase. Under 
P if ip° transforms non-trivially, then P conservation is 
a non-trivial condition that v must satisfy, and its phase 
gets determined. 

Under P we can choose the quarks to transform 
as QiL •<->■ QiR. The Yukawa terms are of the form 



hi.qL 



Qr. 



he with /i,-, real due to CP. Under P if 



we choose (j)° — ?► e'^^cfi" then P invariance of the Yukawa 
term together with real h^^ implies that e^° is real. There- 
fore P = OT (3 — TT for 0° that have Yukawa couplings 
with the quarks. 

Moreover from the P transformation (p" — )■ e'^0° it 
is easy to see that {4>°) = \v\e^^^'^ will conserve P, as is 
needed for VEVs that violate CP. 

Thus the phase of the VEV gets determined in terms 
of the transformation properties of the field under parity. 
For VEVs that have Yukawa coupling with the quarks 
/3 = or TT and so the phase of {(j)°) which is /3/2 becomes 
or tt/2. (3 = corresponds to the case where both P 
and CP are conserved by the VEV and /3 = 7r/2 violates 
CP maximally while conserving P. Both types of fields 
will be present. 

There are a few points we make before we can say 
that the CP violation is also maximal in the quark mass 
matrix in symmetry basis: 

• The phase generated in the quark mass matrix can 
be rotated away if there are only the usual 3 fam- 
ilies. So there is no CKM phase generated unless 
the quark content is extended. This is resolved by 
having a full vectorlike quark family in the model 
as proposed in reference [S* 

• Since there are Higgs fields with both purely imag- 
inary (/3 = tt/2) and real {(3 = 0) VEVs, if they 
contribute to the same element of the quark mass 
matrix then the phase generated in the mass ma- 
trix will not be tt/2 owing to their mixing. We will 
see in the next section that this does not happen 
due to the same symmetries (we call these Strong 
CP Solution Symmetries or SCP for short) that 
are imposed to help solve the strong CP problem. 
So each term in the mass matrix is purely real or 
purely imaginary. Thus a maximal CP phase tt/2 
gets generated in the quark mass matrix in the sym- 
metry basis. On diagonalization to the physical 
basis the real and imaginary parts mix providing a 
CKM phase consistent with the observations. 

• CP can also be broken spontaneously using a CP 
odd, P even real singlet. This case was already dis- 
cussed in reference jS] and its physics is equivalent 
to CP being violated softly by dimension 3 mass 
terms (and with no singlet) as discussed in that 
work. The attractive feature of this case (with or 
without the real singlet) is that no other symme- 
tries (ie SCPs) other than parity and CP need to 
be introduced to solve the strong CP problem. The 
flip side is that it is therefore possible to generate an 
arbitrary (ie non-maximal) CP phase in the mass 
matrix. 

We will estimate lower bounds for and the neutron 
EDM that is radiatively generated in the quark mass ma- 
trix for both the above cases of CP violation.... where it 
is maximal in the quark mass matrix in symmetry basis 



(see section IV A I as well as where it could take on any 
value (see IV B). 



In the next section we set up the P and CP symmetric 
strong CP solving model and show that CP violation is 
maximal in the quark mass matrix. Before we proceed 
a note on the conventions used which are slightly dif- 
ferent from those in this section. So far we defined CP 
transformations so that real VEVs conserve CP. How- 
ever in left-right symmetric theories the convention is to 
define P and CP such that real VEVs of doublet </> con- 
serve P while they could violate CP. We revert to this 
standard convention of left-right symmetric theories in 
what follows, however since the result of maximal CP 
violation are convention independent the particular con- 
vention chosen does not make a difference to the physics. 



III. LEFT-RIGHT MATTER-ANTIMATTER 
SYMMETRIC STRONG CP MODEL 

SU{2)l X SU{2)r X U{1)b-l x P x CP x SCP 



We begin with the left-right symmetric model [T3] 
hasedonGLR = SU{3)cxSU{2)LxSU{2)iixU{l)B-LX 
P with strong CP solving particle content discussed in 
reference [6] where CP was imposed on dimension 4 
terms and broken softly by dimension 3 mass terms, and 
increase the symmetry to include CP as a good symmetry 
of the Lagrangian. Therefore a second Higgs bidoublet is 
added so that it can receive a CP breaking VEV. 

The Higgs fields thus consists of the usual SU{2)ii 
triplet, SU{2)l singlet A^j (and parity related A^) and 
two bi-doublets 4>s and ^a that transform as (1,2,2,0) 
under Glr. The Higgs bi-doublets and triplets are rep- 
resented as complex 2x2 matrices (with the triplet matri- 
ces being traceless) in the relevant iso-spaces, as is usual 
in left-right symmetric model. 

The quark content of the model is the same as in [5] 
where in addition to the usual 3 light chiral families we 
have a vectorlike SU{2)l doublet family Q4L and Q'^ 
that transforms as (3,2,1,1/3), and a parity related vec- 
torlike SU{2)l singlet {SU{2)ii doublet) family Q^r and 
Q'j^ that transforms as (3,1,2,1/3). The prime on the 
quarks is being used in our notation to clearly label the 
mirror component of the vectorlike family. 

Thus there are two heavy top and two heavy bottom 
quarks in the model, one pair from the vectorlike doublet 
and the other from the vectorlike singlet family. Note 
that instead of A^.a, SU{2)l^ii doublets xl,r can also 
be used to break parity and the results in this work will 
apply to them them as well. There are the usual 3 gen- 
erations of leptons and vectorlike lepton families need or 
need not exist. We do not explicitly mention the lep- 
ton sector in this work, other than briefly in the second 
comment of section IVII 



A. P and CP transformations 

Under Parity 5C/(2)i ^ SU{2)r, ^a,s ^ 0l,., Q^L O 
QiR, Q'l ^ Q'r^ ^l <->■ ^R, Note that with these as- 
signments real VEV's for neutral (i.e. diagonal) compo- 
nents of bidoublets (j>a,s will not contribute to P violation. 
QiL,R with i = 1 to 4 correspond to the usual 3 genera- 
tions of quarks and the normal chiral components QiL,R 
of the vectorlike quarks. 

We now choose the CP transformation such that a 
real VEV for (pa will break CP. Under CP we require 



-^ 



QiLM ^ CQ 



iL,Ri 



'L.R 



-^ 



With these transformations note that CP invariance of 
the Lagrangian implies coefficients of all terms (includ- 
ing the Yukawa potential) with an odd (even) number 
of (j)a will be purely imaginary (real). Additionally, the 
Yukawa matrices must be Hermitian (due to P invari- 
ance), dimension 3 quark mass matrix with direct mass 
terms involving vector-like quarks must be real symmet- 
ric (due to CP and P), and dimension 2 Higgs mass pa- 
rameters must be real (due to P). 

We will now discuss a further symmetry that needs to 
be imposed to solve the strong CP problem and see how 
this also leads to a maximal CP phase in the quark mass 
matrix in the symmetry basis. For this we discuss the 
Higgs potential and show that symmetries can be broken 
through real VEVs in the model. 



B. Z2 Symmetry, Higgs Potential and 
P Conserving, CP Violating VEV 

If all parameters of the Higgs potential are real then all 
VEVs that minimize the Higgs potential can be naturally 
real. Note that a real VEV for <j)a will break CP and not 
P and is of interest to us based on arguments of section|IT] 

However there are P and CP invariant quartic terms 
with an odd number of (j)a that have purely imaginary 
couplings and are dangerous for example, 

zrr(A>t0^ + \''cPaUs)Tr{A^j,AR - A[Al) + h.c. (1) 

with A' and A" real and where ipa^s = T-2<?!'a,s'^2- We thus 
need to introduce an additional symmetry under which 
(pa and (ps transform differently such as Z2 with (pa odd 
(or anti-symmetric) and ps even (or symmetric) under 
Z2 to prevent this term. 

Under Z2 we have (pa — > —(pa, QiL — ^ —QiL and 
QiR —>■ —QiR while all other quark and Higgs fields in- 
cluding (ps are invariant under Z2. Note that one gener- 
ation of quarks is chosen to be odd under Z2 to permit 
Yukawa couplings with pa and the remaining generations. 
Z2 ensures that all terms in the Higgs potential are real 
since non-real P and CP invariant terms such as (IT]) van- 
ish due to it. 

Note that any symmetry that sets the purely imagi- 
nary quartic couplings to zero will work as well and we 



designate symmetries such as Z2 as Strong CP solution 
helping symmetries or SCP symmetries for short. 

We will now discuss the Higgs potential in further de- 
tail and show how the symmetry breaking can happen 
to give rise to real Higgs VEVs and a Higgs mass spec- 
trum where all the additional Higgses beyond the stan- 
dard model Higgs are naturally very heavy. 

The most general Higgs potential invariant under 
SU{2)l X SU{2)r X U{1)b-l x P x CP x Z2 is oi the 
form 

Vi„,,=Vt„iAL,AR)+VM{Ps,Pa)+ V V{P„Al,Ar) 



z—a,s 



(2) 
where due to Z2 and gauge invariance there are no terms 
such as in (fl]) that involve both ps and (pa as well 
as A]^ fl, and Vbi contains only quartic terms that in- 
volve an even number of both 0^ and (pa- In the above 
F(02, Al, Aii) + Vtri{AL,AB), with z = a or s, is the 
most general Higgs potential of the minimal left-right 
symmetric model involving only one bi-doublet (p^ which 
is given in references [HI [15] . V, Vtri and Vu are also 
functions of the mass parameters and other coupling con- 
stants of the model and these indices have been sup- 
pressed. All these parameters are real due to CP and 
Z2 that have been imposed. 

We allow for the Higgs potential to have dimen- 
sion 2 terms that break CP and Z2 softly such as 
—^\asTr{p\(ps) + h.c. All the parameters of soft CP and 
Z2 breaking terms such as ^ilas are real due to P. We 
thus have for the Higgs potential 



V, 



Higgs 



Vinv. + ^soft 



(3) 



where the detailed form of Vinv and Vsoft is given in 
Appendix [A] and has all real parameters. 

The neutral component of the Higgs fields can pick up 
VEVs and these are indicated by 



(AL,fl) 




Sl.r 











(4) 



We break SU{2)r x U{1)b-l to C/(l)y in the usual 
way with the neutral component of A^ picking a large 
VEV vji = Sr that sets the scale of parity breaking Mr. 
Due to SU(2)r invariance of the potential vr, can be 
chosen to be real and positive without any loss of gen- 
erality. Note in equation (Al) that the coupling term 



P3Tr{A'fiAR)Tr{A]^AL) with p^ positive has its lowest 
value ii 6l — and therefore this term can ensure that 
Vr, » \6l\. This is the usual way parity breaks in the 
left-right model. 

We now look at the bi-doublet VEVs. Without loss 
of generality we can use SU{2)l invariance to choose n^ 
to be real and positive. The remaining VEVs n's^Ka, k^ 
and 8l can in general be complex and their phases are 
determined by the minimization condition dVuiggs/d't' — 
for all the scalar field components p in the theory. 



We now show that for a scalar potential with all real 
parameters the minimization condition is satisfied when 
the VEVs of the fields are all real. 

We denote with subscripts Re and / the real and imag- 
inary components of every field - for example Ka — 
KaBe + 'ii^ai stc. Vuiggs Can now be written as a sum of 
products of real and imaginary parts of the fields. How- 
ever Vuiggs itself is a real-valued function as the Hermi- 
tian conjugate of every term is also present in it. Con- 
sider a term such as ifJ-ias'^aiK-sRe with just one imagi- 
nary component of the field that can potentially be ob- 
tained from the term —fj,1^gTr{(p\(j)s)- Since /i^^^ is real, 
the Hermitian conjugate —ifJ^ias^^aiKsRe obtained from 
~fil^gTr{(j)l(j)a) provides a canceling contribution. Hence 
there can be no terms with an odd number of the imag- 
inary components of the fields in the real-valued scalar 
potential that has all real parameters. 

This implies that 



dV, 



Higgs 



d(t>. 



•il 



= 



(5) 



01/= 



This is because every term that depends on the imaginary 
components must have at least two of them and when 
the partial derivative is taken and evaluated at the point 
where all the imaginary components are set to zero it 
vanishes. 

Thus we have shown that real VEVs satisfy the min- 
imization conditions. We now show that there exists a 
region of parameter space where the scalar potential is 
at the minimum (and not only at an extremum) for real 
VEVs. 

There are several terms in the Higgs potential that 
depend on the phases of the fields. We choose 
the sign of the parameter of every such term so 
that each of these terms makes a negative contri- 
bution to the Higgs potential when all the field 
values are real and positive. For example we 

choose mass parameters mL> mL> M?as> ^Ls.MLa > 



in (A2) and (A4), and coupling parameters such as 
A2s, A2a, A4s, A4a, has,a2a, a2s, P4 < in (|A2|) and (TASI), 



as well as the ^'s < in (A2) 



With this choice of parameter space it is easy to see 
that the minimum of the Higgs potential will be attained 
when all the fields Ka,s,K'a,B ^'^^ ^l,r pick up real positive 
VEVs (and not any other phases) since each and every 
term of the potential is independently minimized with 
this choice. 

Thus we have proved without any fine-tuning that 
there exists a region of parameter space where the VEVs 
are real. In reality the parameter space of real VEVs 
is much bigger - we need to evaluate the second deriva- 
tives of the scalar potential at the extremum point with 
real VEVs and diagonalize it. The condition that every 
mass eigenvalue so obtained is positive (unless it is zero 
because it is an would-be Goldstone mode) is the only 
condition that the parameters would need to satisfy for 
the extremum to be a minimum. In fact with only one 
bi-doublet it is known that a scalar potential with real 



parameters will only have real VEVs (unless there is a 
lot of fine-tuning) [M] . With two bi-doublets as we now 
have, there maybe some space for non-real VEVs, but a 
significant region of parameter space will give rise to real 
VEVs. The regions with non-real VEVS are phenomeno- 
logically ruled out as they would give rise to a tree-level 
6 and we are left with regions of parameter space with 
real VEVs. 

The remaining minimization conditions 

dVHiggs/d(f>iiie = have to be solved simultaneously to 
obtain the real VEVs in terms of the parameters of the 
potential. 

Note that the natural value for mass parameters (/i) 
of all the dimension 2 terms of the Higgs potential in 
Appendix 1^ is the parity breaking scale Mfj, which is 
heavier than the electro-weak scale and can even be as 
large as the Planck scale. Thus the natural values for 
all the VEVs will be either zero or the scale Mr. To 
get a VEV at the electro-weak scale we need to fine- 
tune parameters in one combination of the minimization 
equations (say in dVniggs/dKs = 0) so that order Mr 
terms cancel giving rise to the weak scale. This is the 
usual fine-tuning that is needed in the minimal left-right 
symmetric model or in any non-supersymmetric theory 
owing to the problem of quadratic divergence. It ensures 
an electro-weak scale VEV for Kg and one light Higgs 
boson which is identified with the standard model Higgs 
boson. 

If we do not allow for any other fine-tuning the other 
Higgs bi-doublet ^^ cannot pick up a VEV in the absence 
of soft CP and Z2 breaking terms in Vgoft ■ 

Therefore we have introduced Vsoft as a perturbation 
so that terms such as n\g^gTr{(l)\4>a) + he will induce a 
VEV to 4>a from the VEVs of (pg . Once we substitute for 
6r in the Higgs potential of Appendix |X] then the terms 
involving Ka for example can be written in the form 



V 



2 ,22 



(6) 



where in order to illustrate the process we have set 
A'ias; /^2sa i^ "^soft to zero. Vbi can be neglected, since for 
electroweak scale VEVs, Vm is much weaker than Vsoft- 
We have also neglected terms involving k^ such as Ka^'a 
for purposes of illustration. In the above miff depends 
on terms such as /xf ^ in the Higgs potential of Appendix [A| 
and the VEV Sr. Further, fJ'ias^''^eff^'^a and Kg are all 
real. 

Minimizing the above potential for Ka gives Ka — 
(Mias/2"T'e//)^s, where the magnitude of term in the 
brackets is < 1 as it was treated as a perturbation 
and can be as low as ~ 10~^ (since Ka is responsible 
for Yukawa couplings between first generation of quarks 
and the rest of the generations). We can see that even 
if /^lasi'Tie// ~ Mji, Ka is naturally small at the elec- 
troweak scale set by k^. 

In general both Ka and k'^ as well as the rest of the 
terms in Vsoft will be present. In the potential terms 
that are quadratic in the fields Ka and k'^ as in equa- 
tion ([6]) will be the only relevant terms in the leading 



order, and will include cross-terms such as Kan'a- Solving 
minimization equations simulatenously these VEVs can 
be obtained in terms of liner combinations of the first 
bi-doublet VEVs Ks and kI^ without any fine-tuning. 

The case with one bi-doublet is considered in detail 
in reference [13]. We have shown how the second bi- 
doublet can naturally pick up VEVs from the first bi- 
doublet VEVs. 

Since there is no fine-tuning other than what is usual 
for the standard model Higgs mass, the natural scale of 
all the other Higgs boson masses is ~ Mfj. 

6l picks up a VEV of the order v'^/vr where v is the 
weak scale, owing to the /? terms in the potential in equa- 
tion ( A2 1 in the Appendix [A] as is usual in the left-right 
symmetric model. The VEV induced for 6l is also natu- 
rally real as the parameters of the potential are real. 

In the next section we look at the fermion mass ma- 
trices which are obtained from the Yukawa matrices and 
real VEVs. 

In place of Vsoft we can also use a complex singlet 
to spontaneously break CP without need for any dimen- 
sion 2 soft CP breaking terms. This is discussed in sec- 
tion IIIIFI In the next section we show how real VEVs 
lead to the solution of the strong CP problem as well as 
to a maximal CP Phase. 



C. No Strong CP Phase, Maximal CP Phase! 

Since Yukawa matrices are Hermitian due to parity, 
real Higgs VEVs imply that quark mass matrices Mu 
and Md are also Hermitian and 6 = ArgDetMuMd van- 
ishes and the strong CP problem is solved at the tree 
level. This can be explicitly seen below from the Hermi- 
tian form of the mass matrices in equation (|9| obtained 
using ^ and ^ 

To note how the maximal CP phase arises recall that 
under CP (f)a —>■ —4>*a: while the rest of the fields trans- 
form in the usual way as given in section |III A[ CP con- 
servation implies that the Yukawa couplings involving (pa 
have to be purely imaginary while the rest of the Yukawa 
couplings are real. Since the VEVs are real this im- 
plies that a purely imaginary contribution is made to 
the quark mass matrix by Yukawas involving (f>a- 

Moreover, the same quark pairs coupling to 0^ cannot 
also couple to <f>s since (pa and (ps transform differently 
under Z2 (or in general under any SCP). That is if the 
purely imaginary Yukawa coupling iQiL(paQjR is permit- 
ted by Z2, for the same i,j we cannot also have the real 
Yukawa term QiL4>sQjR- This implies that after symme- 
try breaking through real VEVs for (pa and (ps , the quark 
mass matrix will not only be Hermitian but will also have 
entries that are either purely real or purely imaginary. 
Thus in the symmetry basis we see that CP violation is 
maximal - ie the phase is 7r/2. It is interesting that the 
vanishing of tree-level 6 (which motivated the Z2 sym- 
metry) implies a maximal phase for CP violation in the 
quark mass matrix in the symmetry basis. 



Concretely if the first generation quarks are odd under 
Z2 and the rest are even, the most general Yukawa terms 
invariant under P x CP x Z2 take the form: 

y^ i[QlL{hij(pa +hlj<Pa)QjR - QjL[h^jPa +hlj<jia)QlR] 
3=2-A 

+ QiL{hn<l>s+hi4>s)QiR + Q'L{!i(t>l+K4>l)Q'R 

+ Ei,i=2-4 QiL{hij(ps + hj(ps)Q]R + h.C. (7) 



where in the last term h. 



_ij h^i, hij — hji due to P, and 
hiij,hij are real due to CP for any i,j in all terms. The 
underline on the Yukawa terms like /i,-,- has been used to 
signify that these Yukawas are in the symmetry basis. 

Due to vectorlike quarks there are also the direct quark 
mass terms 



Mi Q^lQ'r + M*Q'lQ,r + h.C 



(8) 



where the form is due to parity and Mi = M* are real 
due to CP invariance. As we can see in ^ below Mi ^ 
is needed to ensure that the CP phase generated is not 
trivially rotated away. Mi breaks Z2 softly. 

Using ([7]) and (Is]) after electro weak symmetry break- 
ing the up sector quark mass terms are given by 



ith 



M„ = 



/ 


huv 


ih 12V 


^h\,v 


ih'^^v 


Ml \ 




-ihi^v 


K^^v 


ll2,V 


^aV 


M2 




-ih1,v 


h^,v 


h^,v 


K> 


M3 


\ 


—ih^^v 


h^Lv 


hsiV 


KiV 


Mi 


Ml 


M2 


M3 


Mi 


h"'v ) 



(9) 



where ul is the 1x4 row vector (mil, ■■■, u^l) and from 



7|), /i"jW = h 



ij' 



-hijK'a for j ^ 1, h^ljV = kjUs 



htju's 



for i and j 7^ 1 as well as for z = j = 1, and h"v = h Ks 



h Hi'g. Note that {kq,, 



,} are the diagonal elements of 



{<t)a,s) and are real, and u^ = k^ + n'a + k^ + n'g . Note 
that we are using a convention where v ^ 173GeV so 
that a \/2 has been absorbed in its definition. 

The down quark mass matrix Md is similar to the up 
quark mass matrix with u ^ d and K,a,s "^ ''^'a s i'^ equa- 
tions determining them. We use the notation H^^ for the 
Hermitian upper 4x4 sub-matrix of Af„ in ^ divided 
by the VEV v. Likewise H_d is the upper 4x4 submatrix 
of Md divided by v. The underline once again indicates 
that these matrices are in the symmetry basis. 

Note that the above matrix has three light quark mass 
eigenvalues at the weak scale corresponding to the usual 
3 light chi ral fam ilies and 2 heavy quark mass eigen val- 
ues ^ ± \/E|M7F (up to electro- weak corrections) - thus 
there are two heavy down quarks one each from the vec- 
torlike doublet and the vector like singlet family. Similar 
is the case with the up sector. 



D. Role Played by P, CP, and Z2 

To summarize, 

• P implies that the Yukawa m.atrices are Hermitian. 

• CP implies that the Yukawa parameters are either 
purely imaginary (if they couple to cjja) or purely 
real (if they couple to 0s)- 

• In addition to CP, Zi is needed to ensure that there 
are no non-real parameters in the scalar potential. 
Parameters of dimension 2 soft CP and Z2 breaking 
terms are real due to P. Since all parameters of the 
Higgs potential are real all VEVs can be naturally 
real. 

• The quark mass matrices that depend on the prod- 
uct of Hermitian Yukawas with real VEVs are also 
Hermitian thus ensuring that argdet M^Md = 
and solving the strong CP problem. 

• Real VEVs multiplying the purely imaginary or 
real Yukawa terms give rise to purely imaginary 
or real quark mass matrix terms, simultaneously 
ensuring a maximal CP phase in the quark mass 
matrices in the symmetry basis. Here Z2 ensures 
that the same quark pairs do not couple to both 
(/>a and (/)s, keeping the purely imaginary and real 
terms separate. 

• The vectorlike quarks ensure that the CP phase 
generated cannot be trivially rotated away. 

• Real VEVs conserve either P or CP or both. For 
example real VEVS of A^^^ conserve CP while 
breaking P, real VEVs of (/>a conserve P while 
breaking CP, and real VEVs of 0s conserve both 
P and CP. 



E. Spontaneous CP Phase in Fermion Sector 

An interesting feature of equation ([7| is that we can re- 
define the first generation quark fields Qil.r -^ ~iQiL,R 
so that the phase i gets absorbed so that it disappears 
from the Yukawa terms in ([7| and appears in terms of 
the redefined fields in the dimension 3 mass term as 

iMiiQiLQ'R-Q'LQw)- 

Also since all terms of the Higgs potential are real and 
the Higgs VEVs are also real, the model has a complex 
phase only in the dimension 3 fermion mass term and at 
this stage is almost exactly like the model of |£| where 
CP is broken softly by fermion mass terms. The only 
difference is that it is more predictive - CP phase has 
been determined to be 7r/2 and CP violation is actually 
spontaneous since there is a way of defining CP so that 
the Lagrangian is CP symmetric, as that's how we con- 
structed the model. However at the tree-level the spon- 
taneous phase is generated in the quark mass matrix and 
not in the Higgs terms. 



That spontaneous CP violation in multi- Higgs doublet 
models can appear at the tree level in the fermion sector 
instead of the Higgs sector was noticed in reference [TS] 
in a more academic context and its relevance to physical 
things like the strong CP problem and a maximal CP 
phase was not considered. 



F. Domain Wall Problem 

Spontaneous CP violation at the weak scale (due to 
VEV of (pa) can lead to the domain wall problem [17j . 
We have addressed this automatically since we have di- 
mension 2 soft CP breaking terms Vsoft- Moreover, the 
scale of the mass parameters in Vsoft is naturally of the 
order M^. 

However if we do not want any soft CP breaking, we 
can add complex singlet fields aa that can pick up large 
CP violating real VEVs that conserve P, such that under 
CP (J a — > — cr*, under P aa ^ cr^ and under Z2 cTa ^ 
—aa- As before Z2 ensures that dangerous P and CP 
symmetric terms of the type z(Ta(A|jAjj — AJ^A/,) -I- h.c. 
are absent. Further if we change the CP transforma- 
tion of ipa to be (f>a -^ 4>a then the maximal CP phase is 
generated due to the real VEV of aa from the Yukawa 
term ihaaaiQiLQ'n — Q'^QiR) + h.c. The terms such as 
jjilasTr^i^ft) s + h.c. with ^\^g real, now conserve CP while 
breaking Z2 softly. Thus with such a singlet, CP can be 
violated only spontaneously while Z2 can be broken by 
dimension 2 terms and conserved by dimension 3 and 4 
terms. As before Z2 ensures that the real and purely 
imaginary contributions do not mix in the symmetry ba- 
sis in the quarks mass matrix and the CP violation is 
maximal. If aa picks a large VEV higher than the scale 
of inflation then the domain walls will be washed away 
resolving the problem. 

The case of a real singlet has already been considered 
in "^ and it does not require a Z2 symmetry to be im- 
posed as the above dangerous term is automatically ab- 
sent if aa is real and the strong CP problem is solved 
just by P and CP symmetries without needing additional 
help. 

In what follows we allow real ^\as 7^ and CP is either 
broken softly without a singlet, or equivalently there can 
be the singlet and CP is broken spontaneously without 
any soft CP violation. The discussion in the rest of the 
paper does not depend on the details of this choice. 



IV. RADIATIVE CORRECTIONS AND 
NEUTRON EDM 

We follow the procedure in reference (6j where radiative 
corrections to the quark mass matrix and consequently 
to ^ = ArgDetMuMd have been evaluated at the one 
loop level. 

Since P and CP violating terms will interact with one 
another in loops we expect 9 to be generated radiatively. 



Without loss of generality we can choose a flavour basis 
via a common 4x4 orthogonal transformation O on the 
up and down, such that OfjMj — Si/^M where 5n = for 
i < A and ^44 = 1. In the transformed basis Hu and H^ 
are hermitian matrices with complex phases as they are 
obtained from Hermitian matrices in symmetry basis by 
Hu d = O^ H_^ ^O. Correspondingly the mass matrices 

We can now apply a common unitary rotation in the light 
3x3 sector to diagonalize the 3x3 subspace of Hu and 
rotate away all phases in the 4*'' row (and column) so that 
it is real symmetric. In this basis iJ^ is a 4 x 4 hermitian 
matrix with complex coefficients. Writing explicitly the 
matrix elements ft."' of i?„^d, the up and down mass 
matrices look as follows in this physical basis: 



in ( 9 1 transform as M^ , 



M„ = 



fhl.v 








K,v 


\ 





h-,v 





h^^v 











hhv 


hl^v 





K^v 


KiV 


hhv 


K,v 
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Kv ) 



(10) 



Afrf = 



(hi,v 


hi^v 


hi^v 


hf,v 
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hf*v 
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hi,v 
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hiiV 
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(11) 



where in this basis M^ is real symmetric and Md is Her- 
mitian with complex coefficients h'j^ 






We denote by i5M„. . the radiative corrections to the 



Ath 



row and j*'' column of M„ and evaluate the determinant 
to the lowest order in (5M„ and find as in reference [B] 
the most significant contribution to 



argdet(M„ + 5Mu) — Im 



E 

i—u,c,t 



SMu 



SMuii to one loop was evaluated in reference [6] and 
we outline the steps here leading to the same results. We 
keep the discussion more general so that the same results 
also apply to the case of more than one bi-doublet. Note 
that in reference [6] we were working in a basis where 
the 3x3 sector of M^ was diagonal and so we actually 
evaluated the contributions to SMd^^ (instead of SMu^ ) 
in that work. However the procedure is the same when 
we work in the basis where 3x3 sector of M„ is diagonal. 

To generate an imaginary part to SM^^.. we need to 
look for Feynman diagrams where there are both complex 
phases and parity violation. 

In the Feynman gauge the charged would be goldstone 
boson GJ corresponding to SU{2)]^ breaking makes a 
logarithmically divergent contribution to 6 owing to the 



process in Figure [T] Since it is the goldstone mode of 
SU{2)l transformation its Yukawa couplings are [H] 



^L.H. 



u.^ /d d \ u, 

iR ' 4L jR » iL 



FIG. 1. Divergent part of radiative contributions from 
charged would be Goldstone mode and second Higgs doublet 
that are related by parity cancels 



Cy = UiB.h^jdjLG'l - UiLhfjdjRG^ + h.c. 



(13) 



However, since among the Higgses only (jja.s couple to the 
quarks, and they pick up P symmetric real VEVs, and 
the Yukawa matrices are Hermitian, there cannot be a 
net one-loop contribution to 9 from the Higgses in the 
limit that the bidoublets do not couple to A^ - whose 
VEV is the only source of parity violation. 

Therefore the contribution of Gj to 9 will be ex- 
actly canceled by another massless mode in the Higgs 
bidoublets which we can call H2 , and it corresponds 
to a global S'f/(2)^ rotation U^ of the bi-doublets and 
fermions alone. That is, under this unitary transforma- 



tion. 



i>a,sUR 



bUf 



VSQ.R, and Afl 



^R- 



In the symmetry limit A/j Higgses do not couple to the 
bidoublets other than through Tr((/)|0s)Tr AJjA^, which 
due to the trace does not convey SU{2)ii breaking to the 
bi-doublets. 

Now if we turn on the P, CP and Z2 invariant cou- 
pling terms between the bidoublet and A^jj such as 

a3sTr{<f>lcl)sA''„Aii + (l)s(l)lA[AL) + h.c. of equation (|A2| 



(1^) so that U^ is explicitly broken, then when Aji picks 
up a large VEV, i/^ will acquire a large mass M^+ ~ 
y^ct^Mfi, where Mfi is the parity breaking right-handed 
scale. The cancellation of the divergent part of the con- 



tribution from GJ by iJ^ is unaffected by the mass 
picked up by H2 and a finite contribution to SMu^^ is 
generated due to figure [l] whi ch d epends on the product 
{HdHdHy)ii so that 9 due to (12) becomes 



\SM^ 



167r2 



E 

i=\ to 3 

j=l to 4 



Im 



' h'^ h'^ ?7" 
JJ j4''-4i 



\n{M^+/M) 



where M < M^+ and the logarithmic dependence is be- 
cause the divergence of figure [T] beyond M has been can- 
celed beyond scale M^+ . 



9 



As in reference [B] we also note that if M > M^+, 
then since below the scale M we have the left-right sym- 
metric model, parity protects from being generated 
so that the logarithm above is replaced by the factor 
{M ^+ / MY which vanishes in the limit the parity break- 
ing scale M^+ — > 0. 

Note that the bounds on and dn we are calculating 
are on their magnitude and so the sign of the quantities 
is not relevant for these purposes and we are not keeping 
track of it. Also the bounds are an order of magnitude 
estimate and not exact. For example there can be a quark 
mass dependent coefficient of order 1 that multiplies the 
terms within the summation sign of (14). That is, the 
coefficient for j = 4 in ( 14 1 when the heavy quark is in 



the inner loop of figure [Ijcan be different from those for 
j = 1 — 3 for light quarks whose masses can be taken to 
be zero. However since these coefficients are both order 
1 we have combined the two terms instead of writing 
them separately. This simplification does not affect our 
estimates of the lower bounds. 

Keeping these in mind we write in a compact form 



<5A/„ 



167r2 



E 



Im 



1=1 



to 3 

to 4 



' /i<i /id /,« 



(15) 



where the factor 



ln(Mj^+/M) > 1 
(AV/M)2 



for M < M„+ 
for M > M^+ 



(16) 



The choice of the symbol ^ is to remind us that its a 
logarithmic factor for M < M-^+ and for purposes of cal- 
culating lower bounds in this region of parameter space 
can be chosen to be 0(1). We also note that (151 is in 
the basis where 3x3 submatrix of A/„ is diagonal. It can 
also be written in a basis independent manner as shown 
in reference IHl as 



basis of ([9]) can potentially contribute to 9 we can get the 
lower bound by setting maximum number of these ma- 
trix elements to zero while still having sufficient non-zero 
off-diagonal elements to obtain the correct CKM mixing 
angles. 

Since we need non-trivial participation from vector-like 
quarks for CP violation we need to set Mi, M4 7^ 0. In 
addition we need 3 off diagonal Yukawa couplings /i";' to 
obtain the 3 CKM angles if we set M2 = M3 = 0. At 
least one these Yukawa terms (but not all) must be from 
the fourth row or column, or else the the CP phase can 
be rotated away. 

At this stage once we rotate the heavy quarks away and 
go into the physical basis of equations (Il0|) and ( 11 1, the 



light 3x3 quark sector will have 3 off diagonal elements 
(they could be anywhere in up and down sectors) each 
coming from a different non-zero Yukawa of the symme- 
try basis. So a maximal CP phase 7r/2 is transferred to 
the light 3x3 sector and the CKM angle generated is 
too high at 7 - 7r/2. To get the CKM phase of 68.8°, 
we need one more non-zero Yukawa, so finally we need 
a texture with 4 non-zero off-diagonal Yukawa terms in 
the symmetry basis. 

For all such textures we can evaluate 9 and we find the 
lower bound 9 > I0~^^i. For example, let us take all off- 
diagonal Yukawas as well as h^^ in Af„ of([9| to be zero. 
Let us also take h^^ — h^^ — h^i = so that the only 
non-zero off-diagonal Yukawas in Md are Zii2,Zii3:Zi23j 
and ft,34. Also we set Mi, M4 ^ and M3 = M4 = 0. 

We now go into the physical basis by an orthogonal 
transformation in the 1 — 4 plane with si = sin9i — 
Mi/M, ci 



(18) 



0056*1 and M = ^J Mf + Ml we 


obtain 


/ (?^hl^ -SiZiii \ 




Z122 






/i^3 
\-sihli slhli ) 





|5M„ 



1 

167r2 



^ Im(h%h%hl,{hf,\)i)t (17) 



,Z = 1 to 3 

j.,k=l to 4 



where /13X3 is the inverse of the 3x3 submatrix of i?„. 
The contribution to 9 from radiative corrections to M^ 



can be got from (17) with u -^ d. 

Note that the Yukawa terms L|Jt2A/jLjj that give Ma- 
jorana mass to the right handed neutrinos break 5'f/(2)^ 
and a^ cannot naturally be kept small. Therefore we 
expect M^+ ^ Mji. 



A. Neutron EDM With Maximal CP Violation 



We will now use the formulas in ( 15 ) and ( |17[ ) to eval- 
uate the lower bound for 9 in terms of the CKM param- 
eters. As each element of M„ and Af^ in the symmetry 



H, 



j ichj^2 *Ci/ii3 + Sih^, 



34 



Vo 



L!;22 



L±23 
i*33 



lSi/li3 





lSitli2 





/ 

(19) 
where the • in the lower triangle is a short hand notation 
to indicate that the term is obtained from the correspond- 
ing element in the upper triangle by complex conjugation 
to make the matrix Hermitian. By comparing ( 18 1 and 



( 19 ) with equations (^10^ and ( 11 )wc can get the Yukawa's 



in the physical basis (/i"' ) in terms of those in the sym- 



uA\ 



metry basis {h^ ). Now using (151 we find the contribu- 
tion to 9 from radiative corrections of Af„ is; 



1 



167r2 



Im 



uu ud ud 
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uu 

'ni 



10 



167r2 



iu id id 

Ihi 



(20) 



where we have dropped terms of order s^ and taken Ci ~ 
1. Since the Hght 3 x 3 up sector is diagonal in (18), we 

have S1/134 ^ Vublkss ^^^'^ Z*i3/(5iZ*34) ^ cos'f where from 
the Particle Data Group's global CKM fits [IH] we have 
Vub ~ 3.47 X 10"'^ and 7 ~ 68.8° is approximately the 



CKM CP phase. Substituting in ( 20 ) we get 



1 



IGtt^ mf 



v^f,cosje - 10""'^ 



(21) 



where /133 = nih/mt ^ 1/40 is the bottom quark Yukawa. 

We similarly find for all choices of the placement of the 
non-zero off diagonal elements that the strong CP phase 
generated has the lower bound > 10~^^. 

While making the above statement we note that there 
is a special case with the choice h^^ = in equation M 
for alH = 1 to 4 as it can lead to a vanishing 9 at the 
one-loop level. However this texture does not generate 
the observed CKM phase. With Mi ^ in ^ and ^ 
this can be easily seen by first rotating in the 1 — 4 plane, 
then in the 2 — 4 and finally in the 3 — 4 plane and going 
to the physical basis of pO| and (111. Due to the hi- 



erarchy of the Yukawa couplings and Yukawas involving 
the first generation being the smallest, there is no sig- 
nificant mixing generated by the above rotations in the 
light 3x3 sector between Yukawa terms in the first row 
and those in the second and third rows and therefore the 
CKM phase generated is too small. This texture is thus 
not easily viable. 



An interesting thing about equation ( 20 1 is that the 



Yukawas involving the fourth heavy quarks such as h^^ 
are multiplied by the angle si so that the product leads to 
a Yukawa term of the light first generation. This means 
that /134 itself does not have to be small, it is the product 
S1/134 that is small. 

Radiative corrections to Md can also contribute but 
they turn out to be much smaller for textures with max- 
imal CP such as those in (18) and (19 1. However they 



are relevant if the CP phase is not maximal as we will 
see in the next section. 



B. Neutron EDM Bound for General CP Phase 

If CP is broken softly by dimension 3 quark mass terms 
(ie) by choosing Mi to be complex or spontaneously via a 
real singlet (both these cases are studied in reference [5]) 
then the CP phase in the quark mass matrix in the sym- 
metry basis can be non-maximal and arbitrary. We can 
then choose it to be close to the CKM phase 7. This 
means that with just 3 non-zero off-diagonal Yukawa el- 
ements in the symmetry basis (and M2 = M3 = 0) we 
can obtain the needed CKM angles and the CKM CP 
phase. For su ch tex tures, since there is one less non-zero 
term than in IV A we find that the 6 generated will be 



smaller by about two orders of magnitude and accord- 
ingly the lower bound is reduced. 

To see this note that if the CP phase i in the first row 
of equation ( 19 ) is replaced by e*^ with 7' close to the 



CKM phase 7, then we can set hi^ to zero and still obtain 
the needed CKM angles and 7. However if h^^ — then 
the contribution to 9 from equation( 20 1 vanishes. 

On the other hand we find that the contribution to 9 
from radiative corrections to Md remain non-vanishing. 
Since the light 3x3 sector of Md is non-diagonal we need 
to use the basis independent equation |17| with u ^ d to 
estimate them. We find 



1 / ud ud l^d l^u l^u 

-*- T / "-12"-23"'34"-14"-ll 

167r2 \ hdhghh 

'^ 167r2 \ 



hdJiMh 



(22) 



where the product of the Yukawas in the denominator 
corresponds to the determinant of the light 3x3 down- 
sector Yukawa matrix, that is hd^s,b = 'm-d,s,b/iT^t- Sub- 
stituting h^2 = Kis/i22i ^23 ~ ^cbhas Siud as before 
shti = Vusbi^ we get 



1 (mlmb\ . 



mdm^ 



10 



-13.5 



i (23) 



There is one more potential source of radiative correc- 
tions at the one-loop level that is relevant. In order to 
calculate lower bounds we have assumed textures with 
zeros in mass matrices which are not necessarily sup- 
ported by any symmetries and therefore small non-zero 
values will be radiatively generated. To avoid fine-tuning, 
the natural tree-level values for these elements must be 
greater than the one-loop corrections. 

Since the bi-doublets will acquire VEVs diag {k, k'} if 
a matrix element hij is non-zero then the corresponding 

element h^j must be at least of the order hijK' /k (with 
k' << k) since both the up and down sector can receive 
contribution from a term such as hQiL4>QjR- The same 
is also true if the up sector has a non-zero element the 
corresponding down element will also get a small VEV 
by the same argument. 

Even if k' is set to zero it will acquire a one loop cor- 
rection order K/{16TT'^){mb/mt) ~ k/6000 due to the di- 
agram in figure [2J where in the figure ht = h^^, hb = /133 
are the top and bottom quark Yukawas. Thus instead 
of zeros we should have a value for every element in the 
up (down) Yukwas that is at least 1/6000*'* of the corre- 
sponding element in the down (up) Yukawas. Since this 
is very small it doesn't increase the lower bounds we have 
evaluated in (23). 



However where this is relevant is for textures that may 
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FIG. 2. Radiative CP invariant contribution to n' 



have an additional zero in the diagonal element such as: 
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(24) 
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(25) 



/ 



The light 3x3 sub-matrices of the above have a predictive 
texture |19j and for e = it is knovifn in the literature 
as the texture corresponding to solution 5 of Ramond, 
Roberts and Ross (RRR). This texture gives the approx- 
imate predictions Kis = y^nd/rris and Vub — y^n^jjmt. 
If we now evaluate using the above matrices we find 
that the contribution from radiative corrections to M^ 
as well as Md vanishes if e = 0. Thus it maybe possible 
that the lower bound may be weaker than (23) for the 
case where CP violation is non-maximal, but we elimi- 
nate this possibility below. 

As we have argued e cannot be taken to be zero but 
must at least be e ^ e^'^ /i^g/GOOO so as to prevent fine- 
tuning. Evaluating the contribution to 9 from radiative 
corrections to Md by using equation (17 1 with u ^ d we 
get for the above e 

IGtt^ 6000 hdhshb 

sin(27') 



960007r 
-3 X IQ- 






(26) 



where we used {sih24)/hs — Vus, hss ~ ht, Z113//133 ~ 
Vub and hd^s,b,t = ^^d,s,&,t/'^t with the central values 
Vus = 0.2253,' Vub = 'o'.d0347, m^ = 5.1MeP", m^ = 
lOOMeV, rub = A.ldGeV and mt = 173G'ey. 

Note that the above quantity came out higher than 
the previous cases because off-diagonal Yukawas in the 
up quark sector that give rise to any CKM mixing an- 
gles are 10 to 40 times higher than if the CKM angles 
resulted from a down quark Yukawa. And therefore 9 is 
comparatively higher if up sector Yukawas are involved 
in generation of CKM mixing angles. 



Thus we find that the lowest value we can get for 
the case of general (non-maximal) CP violation is given 
by ( 23 ) and this serves as the lower bound. 



bound for 9 is 10 ^^ for non-maximal CP violation and 
10^^^ for maximal CP violation. These bounds corre- 
spond to dn > 10~^^ecm and > 10~^^ecTO respectively. 
For Mj^+ < M there is a further suppression by factor 

(Af^+/M)2 

We now see to what extent we can predict the CKM 
phase that can arise from the maximal phase. 



V. CKM PHASE PREDICTION FROM 
TEXTURE AND FLAVOUR 

Since a maximal CP phase is generated in the sym- 
metry basis it is tempting to see if textures and flavour 
symmetries can help predict the CKM phase that is ob- 
tained from the maximal phase. 



A. Using Texture 

It is well-known that setting some of the diagonal ele- 
ments of Yukawa matrix to zero will lead to a predictive 
mass-mixing angle relation, as the mass would then have 
to come from an off-diagonal element. We can use a sim- 
ilar trick to predict the CKM phase from texture alone. 

As an example we take /i"i = and take M2 , M3 = 
and take the four non-zero off diagonal Yukawas to be 
^121 ^23' ^34 and /i24- We follow the same procedure as 
in section HV Al 

Diagonalizing in the 1-4 plane so that the heavy quarks 
decouple (with s — Mi/M as before), we get the 4x4 
Yukawa matrices in the physical basis to be: 



HiL — 



/o 


sWli 








• 


I122 





C/124 








hi. 





vo 


* 









\ 



(27) 



Hd 



( ich{2 sh^4^ \ 



• 


Z^22 


Z^23 


ishi2 


• 


• 


hi. 


Chi, 


u 


• 


* 






(28) 



We can see from the top 2x2 sub-matrices in the above 
two matrices tha t we get the predictive relations md '^ 
VusTTLs (from (28) with /i22^ = "^s and Vus = c/112/^22) i 
and C0S7 ^ {\/'mu/'mc)/Vus (from (|27 [) with K^2'^ ~ "^c 

and rriu = s^Zi24//^22 and from ( |28[ )) . Thus we get 
7 ~ 79° where we have used the central values r?T,„ = 
2.bMeV, rric = 1290Mel/, Vus = 0.225. 
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If we diagonalize the full 3x3 light quark sub-matrix 
and do the calculation more accurately 7 can be a bit 
lowered to 76°. Comparing this with the experimental 
value of 7 = 68.8° ± 3° we see that while the agreement 
is not too great, the texture does come pretty close and 
gives an idea of how the CKM matrix may get generated 
from the maximal phase. Of course once we turn on the 
terms we set to zero we can easily get corrections to the 
texture prediction and an agreement with experiments. 

We evaluate the 6 for this texture and find it to be 



sin 7 shf2h24ms 



16772 



10""^ 



TOc 



(29) 



We will now see how close we can get to predicting 
the CKM phase using flavour symmetries instead of just 
texture. In the next sub-section we motivate a Z^ sym- 
metry and evaluate how the maximal CP phase and Z^ 
symmetry can help understand how the CKM phase may 
be generated. 



B. CKM Phase Prediction Using Z4 Flavour 

If we take seriously the observed fact that the down 
quark mass md is related to the Cabibo angle via rud = 
V^gms, the two conditions that are needed to obtain this 
relation are: 

• Quark Mass matrix is Hermitian. 

• The down quark mass matrix has a near zero in its 
first diagonal entry. 

Our strong CP solution ensures Hermiticity of the 
quark mass matrices and we impose an additional Z4 
symmetry to ensure the second condition. 

Once Z4 is imposed (j)s and 4>s will transform differ- 
ently under it and therefore both of them will not cou- 
ple to the same quark pairs. This means that ignoring 
off-diagonal elements, the quark mass ratios will be the 
same for all generations since every generation would get 
its mass from 0^ and not from (jjg. Therefore if we fix 
K-'gl Us = rnt/mt ~ 1/40 then since rric/iO and m„/40 are 
much smaller than m^ and irid, they would have to get 
their masses from off-diagonal terms. These constraints 
can lead to the predictions of CKM CP phase, and as 
expected the down quark mass. 

Under Z4 we have 

QlL,4L, Q'n. ^ iQlL,4L,iQ'ii', Q2L,3L — > —'iQ2L,3L', 

Q2R,3R -^ -Q2R,3R] (t>s,4'a "> i(t>s, i'Pa (30) 

while the rest of the quark and Higgs fields are invariant. 
The leptons can also transform non-trivially under Z4 
but we do not write these out explicitly. 

The most general Yukawa potential invariant under 
P xCP X Z2X Z4 is: 



E 

i=l to 4 



h,,QiL<t)sQiR+llQ 



L'Ts'^R 



E ['i-hlj{QlL(i>aQjR- QjL(t>aQlR) 
j=2 to 3 

+ h^i{QjL4>sQRA + Qh4>sQjr)] 

+ih^i{QlL(t>aQR4: " QL4:(f>aQlR.) + h.C. (31) 

where without loss of generality we have chosen a basis 
such that /123 — 0- 

Once the diagonal elements of (/)a_s pick up real VEVs 
diagJKa^s, kJj j} the mass matrix Md is given by 

/ /IuK's i/ii2Ka ihisl^a ikul^'a ^1 \ 



M,= 



^!ll2^a 
V Ml 



h.22^s 



Ui2i^s 
M2 



!l33^s 

M. 



LL2i^s 

M4 



M2 
M3 
M4 



(32) 



Mu is given by a similar matrix with K,a,s ^ f^a s- Note 
that since we set k'^/ks — mb/rrit, M„ is almost diagonal 
with negligible off-diagonal elements while the first two 
diagonal elements of Md are small and need contribution 
from off-diagonal terms to provide the right masses for 
down and strange quarks. 

When we make an orthogonal transformation in the 
2 — 4 plane by an angle S2 ^ M2/\/M| -I- M^ to go to 
the physical basis, the strange quark mass picks up an 
additional contribution so that it becomes 



tK^ 



i.22'^s "1" 2S2/I24 
^^+2s2h 



24' 



± 



rat 



+ 2S2/12 



(33) 



where in the last line we have a ± since rac = Z122''' and 
can be positive or negative depending on the sign of /122 • 
To keep matter simple we set A/3 — 0, h^j^ — 
and we rotate along the 1 — 4 plane by an angle si — 
Mxl^jMl -h Ml + Ml- This takes us to the physical ba- 
sis and the light 3x3 sector of Md becomes: 

ihi2'^a + Sih2A^s ikisf^a + Slh^^Kg 

Md3x3 = I * m° S2Z134KS 

•k -k rUh 

(34) 

We will first do an approximate calculation. If hi^ « 
S1/134 then from (34) we get 

(35) 



S1/S2 ~ \Vub/Vcb\ 



and using equation ( 33 1 with (351 we get 

Vub 



Sl!3i24:^' 



24'^s 



2K 



cb 



m,T 



rricmb 



rrit 



(36) 



Substituting the above in ( 34 1 we get for the phase in the 
down quark mass matrix 



7 



'COS 



'COS 



_1 Slll24:l^s 

Vusm° 

Vub 



2VcbVu. 



It 



mcirLb 



mrrit 



75.7° (37) 
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where we have used from the particle data group's 
global fits to CKM matrix [18] Ks = 0.2253, Vcb = 
0.0410, Kh = 0.00347 and m° ^ m^ = lOOMeV, m^ = 
1.290GeV,mb = 4.19GeV and m* = 173GeV. We have 
used the + from =F to get the lower value of 7°. 

The phase 7° in the mass matrix AI^ (in basis with 
Mu diagonal) is related to the CKM phase 7 roughly by 
the relation 7 = 7° — 3.4° and so we find that due to 
Z4 symmetry the maximal phase n/2 results in a CKM 
phase 



7 = 72°±A(7)r 



where 



A(7)l 



mtVub 



(38) 



(39) 



is the contribution due to a small non-zero h^^. 

There is a subtle point in the above equations. In 
principle /i^^g can also contribute to change the first term 
in (38) from 72° to a different number. This is because 



the exact value depends on equation ( 35 ) which will get 
modified once hi^ is non-zero to 



S2/134KS 
Expanding for small h^-^ we get 



Kb 
Vcb 



si/s2 + 0{h^r^) = Vub/Vcb 



(40) 



(41) 



and thus there is no 0{hi^) correction to the first term 
of equation (38) which remains 72°. Equation (38) 



with (39) is thus correct to 0(^113). 



The other Yukawa we set to zero h^^ can be treated 
analogously and like h^^ it does not change the 72° but 
just makes an additive contribution to the lowest order 
when its turned on. 

Thus we see that the maximal phase generated in the 
quark mass matrix is actually changed to 72° on diago- 
nalization to the physical mass basis due to known terms 
when Z4 is imposed, while there can be a small correction 
from this value by unknown Yukawa terms that can be 
set to zero without affecting anything else. The value 72° 
is obtained by inputting central values from experiments 
for quark masses and CKM parameters (other than the 
CKM CP phase) and therefore it can have an experimen- 
tal error of about 3°. This is thus consistent with the 
current experimentally determined value 7 = 68.8° ± 3° 
even if Yukawas like hi^^ were too small to contribute any 
amount. 

There are a couple of points we need to make before we 
proceed. Equations like ( 35 ) that relate CKM param- 



eters to the symmetry basis Yukawas have been derived 
approximately so they can be presented simply but we 
can do a more exact derivation by using CKM matrix 
to O(A^) in Wolfenstein parametrization. The resulting 



equation ( 38 ) remains the same to within a degree when 
we do this. We also had set M3 ~ and when it is turned 
on it will also make a small contribution to equation ( 38 ) 
and so must be kept small. However including M^ can 
also significantly contribute to 6 and M-^ can be bounded 
based on that. 
We now evaluate 



for this case using (15). We find 



the dominant contribution is due to radiative corrections 
to Mu and is 

^'^(i6^)te)^^^-^-^ 



/ cos 7 



K 



2 2 



3.5 X 10" 



(42) 



We have used Z4 symmetry as an example to show 
how the CKM phase and (i„ may be calculated based 
on fiavour symmetries. Other flavour symmetries can be 
similarly explored. Especially interesting would be those 
that are also motivated from the lepton sector considering 
the neutrino mixing. 



VI. FEW COMMENTS 

Before we conclude we make some comments: 

1. In this work we have shown how P and CP sym- 
metry breaking can lead to predictions for the neu- 
trons electric dipole moment. For a large region of 
parameter space of the model these predictions turn 
out to be in the region that will be experimentally 
probed in the next few years. These experiments [5] 
are often considered to be complimentary to the 
LHC and are currently facing various hurdles [20j . 
However since the neutron EDM generated due to 
P and CP breaking does not get diminished even if 
there is no new physics up until the see-saw, GUT 
or Planck scales these experiments can in fact make 
a discovery due to physics at scales that are out 
of reach for LHC or future colliders. The last re- 
maining parameter of the standard model 9 that 
gives rise to d^, as well as neutrino masses and 
mixing, could both be a consequence of very high 
energy physics emerging from the left-right sym- 
metric strong CP solving model. 

2. Once CP is imposed, it will also apply on the lep- 
ton sector. If there are no vectorlike leptons (anal- 
ogous to the vectorlike quarks) then CP phase is 
not generated at the tree level in the lepton sector, 
and there will be no CP phase in the neutrino mix- 
ing parameters. If there are vectorlike leptons then 
the same mechanism of CP generation can result 
in a maximal phase in the symmetry basis of the 
lepton sector (if CP is broken spontaneously by a 
complex singlet or softly by dimension 2 terms) or 
a general CP phase (if it is broken by real singlet 
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or softly by dimension 3 fermion mass terms). If 
there is Grand Unification, then we would expect 
vector like leptons to also be present. 

3. Our model is generalizable to the supersymmetric 
left-right model (SUSYLR). In fact if we want to 
break CP spontaneously in SUSYLR so as to re- 
duce the number of CP phases such as from the 
gaugino sector \2\\ that could generate too large a 
dn, then we can introduce the vector-like quarks as 
we have done with similar consequences. 

4. We have used a complete vectorlike quark fam- 
ily with SU{2)]^ doublet and parity related iso- 
singlet {SU{2)fi doublet) quarks. In fact if only 
the Higgs triplets A^.a that are the natural candi- 
dates for providing large Majorana masses to the 
right-handed neutrinos are used to break SU{2)]i, 
then this is the only choice of vcctorlike quarks that 
works for solving the strong CP problem, since only 
they can have Yukawa couplings with bi-doublets 
and mix with the lighter quarks. So these quarks 
are uniquely picked in some sense. 

5. The vector like quark masses are protected from 
quadratic divergences from the heavy sector 
since in the model without scalar singlets they 
only couple to the standard model Higgs VEVs 
(through Yukawa couplings with the bi-doublets) 
and through dimension 3 fermion mass terms. Thus 
these quarks can be naturally light including being 
at scales accessible by LHC and future colliders as 
pointed in our previous work 'S^ . Of course the di- 
mension 3 mass terms can also be very heavy like 
being at the GUT scale. Their mass is an indepen- 
dent scale and can be light or heavy. 



Moreover we find that the CP phases of Higgs VEVs 
that are generated spontaneously get determined by the 
P transformation properties since the strong CP solving 
vacuum is such that VEVs that violate CP must conserve 
P. This implies that for Higgses that couple to quarks 
and leptons the CP phases of the VEVs are maximal. 

In addition if CP violation happens either softly by 
dimension 2 terms or spontaneously with a complex sin- 
glet we find that the CP phase generated in the quark 
(in general fermion) mass matrix is also maximal in the 
symmetry basis, thereby providing an understanding of 
the largeness of the CKM phase. The same symmetries 
such as Z2 that are needed to ensure that there is no 
tree-level strong CP phase determine the maximality of 
the weak CP phases in the quark (fermion) mass ma- 
trix. While if CP is broken softly by dimension 3 fermion 
mass terms or spontaneously by a real singlet then the 
tree-level strong CP phase is automatically not present 
and no further symmetries like Z2 are needed and as a 
result the CP violation in the quark (fermion) mass ma- 
trix need not also be maximal. In both cases vectorlike 
quarks (leptons) are needed to ensure the CP violation 
in the quark (lepton) matrices cannot trivially be rotated 
away. 

The strong CP phase is radiatively generated at the 
one-loop level due to the SU{2)]^ breaking charged would 
be goldstone bosons and can be estimated from their 
Yukawa couplings with the quarks that include a full vec- 
torlike quark family. P ensures that there is a mode in 
the Higgses of the right handed sector that cancels the di- 
vergences from the goldstone mode and a finite amount of 
6 is generated that depends logarithmically on the mass 
ratio of the heavy higgs Af^+ and the heavy quark mass 
M for the case where M„+ < M. 

The results for the 6 and (i„ generated are: 



VII. CONCLUDING REMARKS 



• For case of maximal CP violation: 



1950's and 60's presented the surprising discoveries 
that nature through weak interactions violates parity (P) 
and matter-antimatter (CP) symmetries. However while 
doing so nature has left us a puzzling situation of an ex- 
tremely small CP violation by the strong interactions (as 
yet undetected) that seems to imply that there is an un- 
known hidden symmetry, along with a large CKM phase 
that seems to be put in by hand in the form of hard CP 
violation in the Yukawa and weak sector. If the unknown 
symmetries are in fact P and CP themselves which are 
spontaneously or softly broken and therefore hidden from 
our view, then we have shown in this work that the neu- 
tron must know about it since its electric dipole moment 
gets generated when these are broken. 

Finding a neutron EDM consistent with the predic- 
tions in this work can provide evidence for P and CP 
symmetric laws of nature that would have been visible 
at high energies such as those at the time of the origin of 
the universe.. 



> 10"" X 



ln(Af^+/M)>l forM<Mfj+ 



{M^+/My 



for M > Mjj+ 



(43) 



• For case of general CP violation: 



> 10" 
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ln(M^+/M)>l for M < Mfj^ 



{M^+/Mf 



for M > M^+ 



(44) 



and dn can be obtained from the relation dn '^ 2 x 
IQ-^^e ecm [2.^. 

If we are fortunate we will be able to observe the neu- 
tron EDM in the ongoing searches [3] that plan to probe 
dn ^ 10~^^ to 10~^^ecTO in the next few years and up to 
10~^^ecm in the coming two decades. A lot depends on 
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whether M < M^+ or if its larger, then by how much? 
We know that the heavy doublet Higgs mass Mfj+ must 
be around the right-handed symmetry (parity) breaking 
scale. Based on the smallness of the neutrino mass this 
would most naturally be at or above the see-saw scale 
of lO^'^GeV if we assume that the neutrino generations 
have similar Yukawa coupling with the standard model 
Higgs like their counterpart up sector quark generations. 
We also know that M being a fermion mass term is chi- 
rally protected and can be naturally small. Therefore 
there are reasonably good chances to observe the electric 
dipole moment of the neutron in the planned experiments 
that could provide evidence for the hidden P and CP 
symmetries in nature. 



Equations (43) and (44) can also be used in the re- 



verse so as to provide bounds for the mass ratio of the 
heavy higgs and heavy quarks based on non-observation 
of dn- For example the experimental bound < 10 



implies from (W3| for the case of maximal CP viola- 
tion that the vectorlike quarks mass must be M > 
e"i°M^+ -- A/^+/20000 so that the logarithm does 
not contribute more than a factor of 10. This together 
with M^+ ^ lO^^GeV would imply that such vectorlike 
quarks would be out of reach of the LHC if CP violation 
is maximal. However some recent work shows that the 
current experimental bound on 9 could be a factor of 5 
weaker [55] based on the methods used to obtain the re- 
lationship between d„ and 9. If this is the case there is 
still a chance for the vectorlike quarks with maximal CP 
violation to be at the mass scales being probed by the 
LHC. Of course for general, non-maximal CP violation 
there is as yet no relevant constraint on the scale of the 
vectorlike quarks masses from neutron EDM experiments 
and they could be within the LHC reach or at very high 
scales. 



10 



Appendix A: Scalar potential With P x CP x Z2 
and 2 Bi-doublets 

We write down the Higgs potential in terms of the P x CP x Z2 invariant part Vinv , and soft CP and Z2 breaking 
part Vsoft that conserves P. We note from equations (pi) and (pi) that 



V, 



Higgs 



V171V. + Vsoft 



with 



Vinv. = Vt„(AL, Ai?) + T4,(0,, 0a) + J2 Vi(j}„AL,AR) 



where VtriiAi^, Aji) + V{<j>z, A^, Aji) is the most general potential with 1 bi-doublet (j)z as given in reference |14) . 

14„(Al, Afi) ={-^llTriARA^J^) + piTr[(A«A^)]2 + p2Tr(AflAfl,)Tr(A|^A]^) + p4rr(ALAi)rr(A|jA^) + R. ^ L} 

+ P3Tr(AiA[)Tr(A;,Al,) (Al) 



r t 



Vicj,,,AL,An) = - filTric^l^,) - ^lUTri^l^,) + h.c] + Ai,[Tr(40,)]2 + \^,Tt{<j>,,4,X)Tt{<j,Uz) + {A2,[rr(0,4)]2 
+ A4.Tr(0,4)Tr(0,4) + /3i,Tr(0,Afl</)t A[) + /32,Tr(0,Afl4A[) -f /33,Tr(0,Afl^t A[) + h.c.} 
+ {ai,rr(0,0t)rr(AflA|j) +R^L} + {a2zTr{4>,4>\)Tr{ARA^j:{) + R^L + h.c.} 

+ a3.[rr(0t0^AflA],) + rr(0,0t a^a[)] (A2) 



Vsoft = -^J\asTr{(j)l(t)s) - tAasTr{i'l<t>s) - lAsaTr{(i>l4'a) + h.c. 



(A3) 
(A4) 



I 

In the above we have used curly brackets so that terms erations such as i? — )• i or h.c. if mentioned would apply 
that can be obtained by replacing Rhy L or those that to every term and they generate the remaining terms, 
can be obtained by Hermitian conjugation need not be 

independently written. Within each curly bracket the op- Note also that all the parameters in the above po- 

tential are real. a2a,ck2s are real due to CP, and the 
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parameters in Vsoft are real due to P. The remaining 
parameters are all real due to P and they are also real 
due to CP. Either of these symmetries can be used to 
make the remaining parameters real. Z2 ensures there 
are no other terms whose parameters could be non-real. 
Together these symmetries ensure that all parameters of 
the Higgs potential are real. 

Vbi is made of several terms that necessarily involve 



both (pa and (ps- -^2 and gauge symmetry implies that 
each bi-doublet will occur twice in Vm- In the above we 
have only written some of the terms of Vti- 
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